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Using multi-dimensional particle-in-cell (PIC) simulations we study ion acceleration from a foil 
irradiated by a circularly polarized laser pulse at lO^^W/cm^ intensity. When the foil is shaped ini- 
tially in the transverse direction to match the laser intensity profile, the center part of the target can 
be uniformly accelerated for a longer time compared to a usual flat target. Target deformation and 
undesirable plasma heating are effectively suppressed. The final energy spectrum of the accelerated 
ion beam is improved dramatically. Collimated GeV quasi-mono-energetic ion beams carrying as 
much as 18% of the laser energy are observed in multi-dimensional simulations. Radiation damping 
effects are also checked in the simulations. 

PACS numbers: 41.75.Jv, 52.38.-r, 52.38.Kd 



Ion acceleration by ultraintense ultrashort laser pulse 
interacting with solid targets has been extensively stud- 
ied in the last decade [li [3, 0, H] due to a number of 
prospective applications, such as proton therapy [5|, pro- 
ton imaging [6], ion beams ignition for laser fusion f?], 
etc. Recently along with the progress of plasma mirror 
technology, the ion acceleration from laser-foil interaction 
has attracted much more attention. It has been shown in 
one-dimensional (ID) particle-in-cell (PIC) simulations 
that specially by use of circularly polarized (CP) laser 
pulses monoenergetic ion beams can be generated in prin- 
ciple (l, H, [13, El- The key effect here is the suppression 
of electron heating [l2|, which otherwise would dis- 
perse the plasma electrons in space and destroy the mo- 
noenergetic acceleration. When a CP pulse is used, the 
laser ponderomotive force does not oscillate and mainly 
pushes the electrons forward. A static electric field is 
formed at the laser front. Ions in the region of electron 
depletion are then accelerated by this static field. The- 
oretical models and simulations based on ID geometry 
have shown a very promising scaling law for the final en- 
ergy spectrum of the accelerated ions. It predicts quasi- 
monoenergetic GeV ion beams for sufficiently long driver 
laser pulses @, d, [l3|- This kind of acceleration belongs 
to the laser pressure dominated acceleration (LPDA). 

However, multi-dimensional simulations show that the 
acceleration structure is not so stable Electrons 
and ions are inevitably dispersed transversely when the 
target is deformed and heated by the driving pulse. Be- 
sides target deformation, instabilities are also another 
fatal problem. The transverse instability of the acceler- 
ating structure limits the maximum energy of the accel- 
erated ions and broadens the final energy spectrum. To 
overcome this, we have recently studied the laser mode 
effects on the acceleration structure and proposed to use 
combined laser pulses for ion acceleration and collima- 



tion In this letter we restudy the problem by con- 
sidering the target shaping, which might be easier to re- 
alize in experiments. The target fabrication has already 
been applied before for the ion acceleration in the target 
normal sheath acceleration (TNSA) regime, where the 
optimization of ion energy spectrum can be achieved by 
target shaping (l5Ul6lJl7| . The experiment with a micro 
structured target [18] produced a mono-energetic proton 
beam in the TNSA regime. 

In the LPDA regime, the target has usually a thickness 
of a few hundred nanometers. Fortunately, due to the 
rapid progress in the nano-technology, structured nano 
thickness targets can be engineered today. In the present 
study, using 2D- and 3D-PIC simulations we show how 
to optimize the collimation and mono-chromaticity of the 
accelerated ion beams via the target shaping. By optimal 
matching, a collimated, GeV quasi-monoenergetic proton 
beam can be generated by a CP laser pulse at lO^^W/cm^ 
intensity incident on a shaped foil target (SFT) with the 
thickness of a few hundreds nanometer. 

Firstly we study the target deformation under the in- 
teraction of a laser pulse. From the momentum conser- 
vation law between the laser pulse and the target, the 
evolution of the target area momentum {p) can be de- 
scribed as following (il. [l3|: 
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where / is the laser intensity, a is the target area density. 
For the velocity evolution of the target one obtains: 
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Where E indicates the intensity of laser electric filed, no 
and lo are the target initial density and thickness, respec- 
tively. It shows that the energy spread of accelerated ions 
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Figure 1: (color online) (a) Layout of shaped target, (b) 
Energy spectrum of ions at t = 30To and t = 40To. Here 
To = 3.33/s represents the laser period. The orange dashed 
line shows the exponential decrease of the spectrum, (c) En- 
ergy evolution of accelerated ions from multi dimensional PIC 
simulations and ID theoretical calculation. Here t irr repre- 
sents the time of laser irradiation on target, (d) Energy spec- 
trum of electrons at t = 40To in the simulations with and 
without the radiation damping effect. 



depends on the transverse variation of the local ratio of 
laser intensity to the target area density. The distance 
the ions pass in the target is: 5(r)(x£^^(t, x, r)/^^. An 
initially flat target is inevitably deformed, if the laser 
intensity is not uniform transversely. The target defor- 
mation quickly destroys the acceleration structure and 
deteriorates the beam quality. From Eq. [2] we see that 
a target can be kept flat if its areal density a is shaped 
properly. For the usual transversely Gaussian pulse, one 
can use a target with the Gaussian thickness distribution 
as shown in Fig. [TJa). In the following simulations, the 
distribution of the target thickness along the transverse 
direction is: 



/ = max{h,loX exp[(-r^/cr|^)'^]} 



(3) 



Here r is the transverse distance to the laser axis, 
hjoi o-T^m are the shape factors, which are shown in 
Fig. Ha). 

First, we use 2D simulations to flnd the optimal pa- 
rameter region because they are computationally less ex- 
pensive than simulations in the full 3D geometry. The 
total simulation box is 32A(x) x 32A(?/), here A is the 
laser wavelength. The foil plasma consists of two species: 
electrons and protons. They are initially located in the 
region 5X < x < 5.3A with the density of n = lOOnc, 
where ric = a;^me/47re^ is the critical density for the 
laser pulse with the frequency uj. For 1/im laser pulse it 
is nc = 1.1 X 10^^/cm^. Here, we present the results for a 
shaped foil target whose parameters are Iq = 0.3A, = 
7A,/i = 0.15A, m = 1. For the flat target, we just set 




Figure 2: (color online) Spatial density distribution of ions in 
the flat target case (a) and in the SFT case (b) at t = 25To. 
Spatial distribution of the laser intensity {Ey-\-Ez) in the flat 
target case (c) and in the SFT case (d) at t = 25To. 



h = 0.3A, other parameters are the same. Thus, the to- 
tal number of ions in the center part of the SFT is orig- 
inally less than that in the case of the flat target. The 
target is shaped along the Y direction in the 2D case and 
in the radial direction in the 3D case. The normalized 
maximum amplitude of the laser electric fleld at the focus 
is ao = eEo/mujc = 100. This corresponds to the laser 
intensity of / = 2.76 x 10^^ VK/cm^ for the assumed wave- 
length A = 1/im. The fufl width half maximum (FWHM) 
radius of the focal spot is (Jl = 8A. The laser pulse has 
a trapezoidal temporal intensity proflle (linear growth - 
plateau - linear decrease) , with lA/c — 8A/c— lA/c. Thus, 
the total laser pulse energy is about 824 J. At t = the 
laser pulse enters the simulation box from the left bound- 
ary. Since the normalized laser electric fleld we used here 
is about 100 and the laser pulse is circularly polarized, 
we include the electron radiation damping in the VLPL 
code and check its effects on the electron cooling (i9I.[2q|. 

Fig. [T^b) shows the energy spectrum of the acceler- 
ated ions at t = 30To and t = 40To for the flat and 
shaped targets in the 2D-PIC simulations. The flat tar- 
get produces no obvious peak structure in the spectrum. 
Instead, the spectrum shows an exponential decrease like 
dN/dEk - 10^ X exp[-^fe/^eff] with E^^ ^ 500 MeV for 
Ek > 300MeV and a cutoff energy 1.7 GeV at 40To. 

When a SFT is used with the transverse shape factor 
(Tt = 7A, the spectrum becomes quasi-monoenergetic. 
The energy of the peak is about 880MeV at t = 30To and 
later reaches 1.2GeV at t = 40To, which are very closed 
to the analytical values obtained by solving Eq. [2l The 
analytical values are shown in Fig. [Tfc) by the solid line. 
As we can see the maximum ion energy at tsim = 40To 
in 2D simulation is a bit higher than the ID theoretical 
value. This is because of reduction of the target area 
density during the interaction. Although the maximum 
cutoff energy of the ions in the SFT case is lower than 
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that in the flat target case, much more protons are ac- 
celerated in a much narrower region, which beneflts the 
further application of the accelerated proton beams. 

To show the polarization effect, a linearly polarized 
laser pulse is used. The magenta dashed line in Fig.IlJb) 
shows the ion energy spectrum at t = 40To. In this case, 
the electrons are easily heated and scattered by the oscil- 
lating part of the laser ponderomotive force. The target 
becomes transparent to the pulse very soon. Ions are only 
accelerated by the spatially dispersed electron cloud and 
cannot get as high energy as in the CP pulse case. The 
spectrum is again exponential with a lower cutoff. 

Figures [2](a,b) show the spatial distribution of ions at 
t = 25To in the two compared simulations. The target 
shaping leads to a more transversely uniform ion accel- 
eration. The initially flat target in contrary is deformed 
and a natural cone builds up during the interaction. The 
laser intensity distribution shown in Fig. [2Kc,d) conflrms 
this. The natural cone focuses the lateral laser energy 
to the center and thus reinforces the on-axis ion accel- 
eration. On one hand, this effects destroys the foil, but 
on the other hand it leads to the higher cutoff energy as 
shown in Fig. W[h). Besides this, when the laser pulse 
irradiates the cone, electrons are easily extracted by the 
laser fleld out from the inner wall of the cone and heated 
because of the oblique incidence. These heated electrons 
disperse in space and pollute the acceleration structure, 
which destroys the mono-energetic character of the ion 
spectrum. By use of the shaped target, these undesir- 
able effects are reduced dramatically. 

The angular distributions of the accelerated ions in 
the two target cases are presented in Fig. [3l It shows 
that in the SFT case the accelerated ions mainly move 
forward. However, in the flat target case, only a small 
portion of the highly energetic ions moves forward. Ions 




Figure 3: (color online) Angular distribution of ions at t = 
30To and t = 40To. (a) and (b) correspond to flat target case; 
(c) and (d) correspond to SFT case. The color represents the 
relative ion number. 



in the middle energy range get a considerable transverse 
momentum. From the simulation, we flnd the average 
emission angle for the ions whose energy is larger than 
1 GeV is about 2.7° in the SFT case and 5.22° in the flat 
target case. The number of ions in this energy range is 
1.9 times larger in the SFT case as compared with the 
flat target. Clearly, both the collimation and the total 
flux of accelerated ions are improved in the SFT case. 

To ensure that these effects are not a 2D artefact, we 
perform full 3D simulations. For the shaped target, we 
use (Jt = 6A in the 3D-simulation. This value appears 
to be the optimal one for the assumed laser parameters. 
The initial position of the target is moved to x = 2A to 
reduce the computational cost. The laser longitudinal 
proflle is also reduced to be :lA/c — 5A/c — lA/c. Other 
parameters are the same as those in the 2D simulation 
above. The electron and ion distributions at t = 20To are 
shown in Fig. [H As we see, in the SFT case a compact 
target sheath breaks out from the rest of the foil and 
is accelerated by the laser pulse. The same effect was 
also observed in the 2D simulation as shown in Fig. [2Kb). 
In contrast, in the flat target case. Fig. SJc) displays a 
continuously dispersing ion density distribution. The ion 
energy spectrum shown in Fig. S^d) also conflrms the 
quasi-monoenergetic peak in the SFT case. The number 
of the ions with energy larger than SOOMeV are 5.09 x 10^^ 
and 6.63 x 10^^ for the flat target and shaped target, 
respectively. And their energies are 5.05 x 10^^ MeV and 
6.16 X 10^^ MeV, the conversion efficiencies are 14.72% 
and 17.94%, respectively. It deserves to note that in the 
3D case the simulation results are far larger than the ID 
analytical values as shown in Fig. HJc). The calculated 
peak value of the ions is 635MeV at tsim = 20To, however, 
the simulation result is 910MeV. This larger difference is 
also due to the target dispersion, which is stronger in 
the 3D case. Generally the ID estimation gives a higher 
energy conversion efficiency and a lower peak energy. 

We also change the parameters h and ar to see their 
effects on the final spectrum. If we fix other parameters 
and vary li but keep li < 0.2 A, the mono-energetic part 
of the final spectrum almost does not change, only the 
lower energy part increases with li. This means that only 
the ions in the center part of the target contribute to the 
final mono-energetic peak. Correspondingly, the target 
width (Jt is a critical parameter for the final spectrum. 
We find in the 2D geometry that when ar/cTL^O.Sb, the 
optimum spectrum appears. For the present simulation, 
when ar/cTL ^ [0.4,1.0] the monoenergetic peak exists. 
In the 3D geometry, the optimum value of ar is usually 
smaller than that in the 2D case. Detailed analysis and 
simulation work on the target shape parameters such as 
m^li^dT is still continuing. 

Further, we perform simulations to check the radiation 
damping effects, which was found to be very important 
for transparent nano targets [2l|. In our simulation pa- 
rameter region (a < 100, no = 100 and /i,/o ^ 0-1), 
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Figure 4: (color online) Spatial density distribution of elec- 
trons (a) and ions (b) in the 3D simulation at t = 20To. (c) 
Spatial density distribution of ions in the flat target case, (d) 
Energy spectrum of ions in the SFT case and flat target case 
at t = 15To and t = 20To. At t = 20To the driven laser pulse 
has already left away from the target. 



In conclusion, by target shaping we have improved 
both the maximum peak energy and collimation of the 
accelerated ions. The shaped target suppresses the tar- 
get deformation, electron heating and makes the ion ac- 
celeration much more uniform in the transverse direc- 
tion as compared with the plain flat target. Because of 
the absence of laser focusing by the natural cone, the 
maximum cutoff ion energy is smaller in the shaped tar- 
get case. However, more ions are concentrated in the 
quasi-monoenergetic peak. Radiation cooling effects are 
checked and they are not so important in the present 
laser-plasma configuration. 
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a = 3.33no/o ^ tttioIo, which satisfies the optimum 
matching of target thickness and laser intensity [22], the 
target is not transparent to the laser pulse and electrons 
are not completely exposed to the laser field. Most of 
the electrons are inside the target with relatively low 
transverse momenta. Our simulations show that only 
the highest energy electrons are cooled down as shown in 
Fig. Hfd). The x-ray photons are mainly radiated at the 
angle of {0 = 30^ — 40^) and the average photon energy 
is about lOkeV. We find that in the present case 21m J 
energy is transferred to the x-ray photons. For the ion 
acceleration, the very high energy electrons are not so 
important, only the bulk of the ponderomotive energy 
electrons constructs the local charge separation field. So 
in the present case the radiation cooling effects are not 
strong enough to benefit ion acceleration. Our simulation 
results also demonstrate there is no observable improving 
on the final ion spectrum when the radiation damping is 
considered. 

In addition, we should point out that the target shap- 
ing only helps to suppress the target deformation and 
electron heating. The surface instabilities still exist. As 
we can see from Fig. [2Kb) and Fig.SJc), plasma bunches 
have been formed in the plasma front due to these in- 
stabilities. They make a little peak in the ion spectrum 
in the earlier periods. However, they evolve further if 
the laser pulse is long enough and make electron heating. 
Suppression of such kinds of instabilities should be an 
important work both for the laser ion acceleration itself 
and for the fast ignition of inertial fusion targets based 
on laser-accelerated ion beams [7]. 
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